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Abstract 

We revisit the problem of the critical velocity of a clean one-dimensional superconductor. At 
the level of mean-field theory, we find that the zero-temperature value of the critical velocity — 
the uniform velocity of the superfluid condensate at which the superconducting state becomes 
unstable — is a factor of ^/2 smaller than the Landau critical velocity. This is in contrast to a prior 
finding, which held that the critical velocity is equal to the Landau critical velocity. The smaller 
value of the critical velocity, which our analysis yields, is the result of a pre-emptive Clogston- 
Chandrasekhar-like discontinuous phase transition, and is an analog of the threshold value of the 
uniform exchange-field of a superconductor, previously investigated by Sarma and by Maki and 
Tsuneto. We also consider the impact of nonzero temperature, study critical currents, and examine 
metastability and its limits in the temperature versus flow-velocity phase diagram. In addition, we 
comment on the effects of electron scattering by impurities. 
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I. INTRODUCTION 

The Landau criterion [ij] concerns the threshold velocity f l of an obstacle moving through 
stationary superfiuids at zero temperature, beyond which excitations are created and super- 
fluidity is lost. By Galilean invariance, the criterion also applies to the threshold velocity 
vl of a uniformly flowing superfluid; see, e.g., Ref. Q]. The value of the threshold velocity 
vl is determined by the following formula 

Vl = min f — ) , (1) 

where Ep is the excitation spectrum. For superconductors the Landau criterion gives 
Vl = Aq/ pp, where 2Ao is the pairing gap at zero temperature and zero flow, and pp = hkp 
is the Fermi momentum of the entities that are paired (with kp being the associated angular 
wavenumber). However, it was found by Rogers that in three spatial dimensions 

the critical velocity v^ — the uniform velocity of the superfluid condensate at which the su- 
perconducting state becomes unstable — of a clean superconductor exceeds vl- Furthermore, 
as discussed, e.g., in Ref jsl, for superflow velocities v in the range v^ < v < Vc, gapless 
excitations occur in clean superconductors. Moreover, the ratio v^^^/vl has been found 
to be e/2 1.359). In contrast, it was found [sl that in two dimensions v^^^ /vl = 1. 
As for the case of one dimension, Bagwell reported that vi^^^ /vl = 1, as found in two 
dimensions. In contrast to the case of three dimensions, in neither one nor two dimensions 
is gapless superconductivity predicted to occur in the presence of flow. 

In this Paper, we re-analyze the critical velocity of a clean superconductor in one spatial 
dimension via a mean-fleld treatment, and obtain the dependence of this velocity on the 
temperature T, along with the associated temperature- velocity phase diagram. In particular, 
we flnd that at zero temperature the critical velocity is smaller by a factor of a/2 than 
the Landau critical velocity vl. Even though for v = vl/V2 a gap in the quasiparticle 
excitation spectrum remains, the superconducting state flrst becomes unstable there. This 
is in contrast to a previous report [6|, which held that the critical velocity is vl, i.e., the 
velocity at which the gap in the quasiparticle spectrum closes and excitations proliferate 
so as to destroy superconductivity. The smaller value of the critical velocity, which we 
obtain here, is the result of a pre-emptive Clogston-Chandrasekhar-like discontinuous phase 



transition [7| , and is analogous to t 
previously investigated by Sarma 



le threshold uniform exchange- 
and by Maki and Tsuneto 



ield in a superconductor. 



At low temperatures 



(i.e., for T < T* 0.556 T^f, where is the critical temperature in the absence of flow), 
the transition from the superconducting to the normal state, occurring due to the presence 
of flow, remains discontinuous. In contrast, for T > T* the transition is continuous, i.e., 
it is associated with a continuously vanishing order parameter, e.g., as T approaches the 
flow- velocity-dependent transition temperature Tc{v) from below. 

From the experimental standpoint, it is usually more convenient to control the condensate 
current rather than its velocity. (An exception is provided by a closed loop threaded by 
magnetic flux.) The more widely appropriate physical observable is thus the critical current, 
i.e., the maximum equilibrium current that a superconductor can sustain, and this is what 
experiments on superconductors frequently measure. It should, however, be remarked that 
recent experiments on t rapp ed condensates of atomic gases make it possible to probe the 



io|. 



critical velocity direct 

Langer and Fisher [llj introduced a fresh perspective on the issue of the critical velocity 
of a sueprfluid. They explained that — as a matter of principle — superflow is inherently un- 
stable, even at velocities below those that allow quasiparticle excitations to proliferate. This 
instability results from the possible occurrence of intrinsic, topologically allowed, collective 
excitation events, in which, e.g., a vortex ring nucleates locally, traverses an Arrhenius en- 
ergy barrier, and grows without bound, thus eradicating a quantum of flow velocity. In 
the case of narrow channels, such as the one- dimensional superconductors discussed here, 
the fluctuations take the form of phase slip events, which can be either thermal or 



quantum 



131 ]. But, before one takes into account the effects of fluctuations, it is important 



to understand and clarify the behavior at the mean-field level, which is mainly what we 
consider in the present Paper. In other words, the work reported here assumes that the 
rate of such toplogical fluctuations is negligibly small, so that our critical-velocity results, 
established on the basis of thermodynamic stability and the possibility of quasiparticle pro- 
duction, retain a use. In the present Paper we focus on issues of phases of thermodynamic 
equilibrium, competitions between them, and metastability; we do not attempt to address 
issues of kinetics, such as the rates at which phase transitions proceed and metastable states 
evolve into stable ones. 

The remainder of the Paper is organized as follows. In Sec. [Ill we diagonalize the Hamilto- 
nian of a superconductor in the presence of flow, and in Sec. Illll we calculate the free energy 
of this system. In Sec. II VI we use the self-consistency of the order parameter to calculate the 
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dependence of the order parameter on the flow velocity at various temperatures. In Sec. |V] 
we calculate the dependence of the transition temperature on the flow velocity, assuming 
the transition to be continuous, i.e., in the limit of linear instability. Double solutions for 
the transition temperature turn out to exist for larger values of the flow velocity. We resolve 
this issue in Sec. I VI I by identifying the globally stable solution, i.e., the one that corresponds 
to the lower free energy. In Sec. IVIII we address the existence of metastable solutions, and 
discuss the extent to which these provide a superflow-based analog of the phenomena of su- 
percooling and superheating. We also obtain the analog of the superheating limit, together 
with the equilibrium phase boundary and the supercooling limit. Having ascertained the 
true, equilibrium, velocity-dependent order parameter, we calculate in Sec. IVIIII the depen- 
dence of the supercurrent on the flow velocity, and hence determine the critical current and 
superfluid density. In Sec. IIXI we briefly discuss the effect of disorder on the nature of the 
transition. In Sec. |X] we contrast the critical- velocity results for one dimension with those 
for two and three dimensions. We conclude in Sec. IXIi 



II. HAMILTONIAN 



We shall use a microscopic approach to discuss an effectively one-dimensional super- 
conducting system, in which there is spatially uniform flow at velocity v = q/m. Here, 
m is the electron mass and q is the corresponding momentum associated with the flow 
(or equivalently, the angular wavenumber, as we shall set h = 1). (What we mean by 
"effectively one-dimensional" is that the transverse dimensions are much smaller than the 
zero-temperature coherence length fo-) We begin by writing down a Hamiltonian that is 
equivalent to Eq. (HI. 4) of Rogers 3] and Eq. (1) of Nozieres and Schmitt-Rink 

^ j 4+9,(7 CA.-+g,<7 + ^ Vfe,,fc2 4i+g,T ^-ki+q,l ^-^2+9,^ Cfca+g,! > (2) 

where Cka and c|.^ are respectively the annihilation and creation operators of electrons of 
spin-projection a =t or | at (one-dimensional) momentum k, Vuk' is the BCS pairing 
potential and /i is the chemical potential. Following Rogers 3], we assume that Vkk' 
depends on the difference between k and k', and hence is independent of q. To justify our 
ignoring any q dependence, we note that the critical momentum q is no bigger than roughly 
mvL, and this is roughly {a/^o)kF, i-e., several orders of magnitude smaller than kp, where 



a is the atomic spacing and is the zero-temperature superconducting coherence length. 
This form of Hamiltonian anticipates that any Cooper pairs are formed via the pairing of 
electrons in states {k + q, |) and (— /c + g, J,), and that the resulting pairs have center-of-mass 
momentum 2q. 



Next, we make the mean-field approximation and adopt the BCS form for Vkk' 15j, thus 
arriving at the Hamiltonian 

^ = ^ + 4+g,T Cfc+<?,T + ^q{-k) cLfc+g,| C-k+q,l 



in which g = \Vkk'\ is the magnitude of Vkk' in the momentum range for which the pairing 
potential is nonzero, eg{±k) = (g ± kY /2m — /i, and A = — Vk^k'{c_k'+q,i Cfc'+g,|) is the 
self-consistency condition on the order parameter (and similarly for A*). This mean-field 
Hamiltonian can be solved via the Bogoliubov-Valatin transformation 

Ck+qA = Uk'Jl-k + Vkltk and C-k+q,l = Mfc72;fc - Vk^lk , (4) 



with 

where, for convenience, we have made the definitions Sk = {£q{k) + eq{—k)) /2 = k'^/2m — 
{fi — q^/2'rrij and Ek = ^/^l + l^p. We choose A to be real and non-negative, and thus may 
drop the absolute value on A in the definition of Ek. With this procedure, the Hamiltonian 
becomes diagonal: 

H = — + J2 {{Ek + kv)^l,^,.,k + {Ek - kv)^l,^2;k) + - Ek). (6) 

^ k k 

This expression also holds for higher dimensions, provided we replace k hy k, v hy v., and 
kv hy k ■ V . 

III. FREE ENERGY 

As the Hamiltonian has been diagonalized, we can readily evaluate the Helmholtz free 
energy F of the system: 

F={H)-TS={H)+TY, [/l;fcln/l;fc + (l-/l;fc)ln(l-/l;fc) + /2;fcln/2;fc + (l-/2;fc)ln(l-/2;fc) 

k 

(7) 



where /i and /2 are Fermi distribution functions, defined via 
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(■ ■ ■ ) indicates an average, weigfited by the equihbrium density matrix, and /3 = 1/T is the 
inverse temperature (with Boltzmann's constant set to unity). We therefore arrive at the 
result: 



One may equally well calculate the partition function Z, so as to obtain the free energy via 
F = —kBTln{Z). The order parameter is to be determined self-consistently, via 



We note that because the momentum sum runs over both positive and negative values one, 
can safely replace fi-k by f2-k (or vice versa) in this equation, and thus arrives at the result: 



which is identical to Eq. (21) of Ref. [6|. This equation can also be derived by demanding 
stationarity of the free energy, i.e., 6F/6A = 0. The above results, Eqs. ([7j) to (ITT!) hold for 
two- and three-dimensional systems as well, provided we replace k hj k, v by v, and kv by 
k ■ V, but in the following we shall focus mainly on the case of one dimension. 

IV. ORDER PARAMETER 

From Eq. (fTOll . or equivalently Eq. (fTTl) . we can determine the values of the order pa- 
rameter at any temperature and superflow velocity. Of course, there is always the trivial 
solution A = 0, which corresponds to the normal state. However, in this section our focus 
is on nontrivial solutions. 

A. Zero temperature 

Let us first examine the limit of zero temperature. In this limit, the self-consistency 
condition ( fTTl) becomes 






(10) 




(11) 




(12) 
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FIG. 1: Self-consistent solutions for the order parameter A (normalized to the zero-temperature, 
zero-flow value Aq) vs. superfluid velocity v (in units of vl) at various temperatures T = 
(0,0.223, 0.445, 0.668, 0.890, 0.980) r° (from top-right to bottom-left). Note the multivaluedness 
of A, which occurs for the lowest three values of the temperatures. It turns out that the lower 
branches of the order parameter have the maximum in the free energy compared to the upper 
branches and the trivial A = solution. 

where 6(x) is the Heaviside step function. By linearizing the fc-dependent spectrum around 
k = kp, and exchanging g for the (T, v) = (0, 0) value of the order parameter Aq, we arrive 
at the following result for A(T, i;)!^^^: 

'0 



(13) 



This condition yields two branches of solutions for the order parameter: (1) If kpv < A, we 
obtain A = Aq; on the other hand, (2) if kpv > A, the condition becomes 



In 



7? 



A2 



sinh 



-1 Vk 



2„2 



A2 



(14) 



which leads to A^ = 2A;^ v Aq — Aq. These two zero-temperature solutions were first obtained 
by Sarma 8|] in connection with the exchange-field effect in superconductors. In Fig. [1] 
we show these two zero-temperature solutions [i.e., the horizontal and parabolic curves 
emanating from the upper right point (1, 1)], together with the corresponding solutions at 
nonzero temperatures. By continuity, it is not surprising that there is a low-temperature 
regime in which there is multivaluedness in the solutions for A, as we shall see in the following 
subsection and in Fig. [H 
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FIG. 2: Critical temperature (from linear instability) Tc{v) (in unit of T^) vs. superfluid velocity v 
(in units of vl), as obtained from Eq. (I17p . Note the occurrence of the unphysical multivalueness. 



B. Nonzero temperatures 

For these, we can solve Eq. ffTOj) . or equivalently Eq. ffTTj) . numerically, and obtain the 
velocity-dependent order parameter A{T,v) at arbitrary temperatures, as illustrated in 
Fig. [1] for several values of the temperature. We remark that in solving for the order param- 
eter we have linearized the spectrum Ek ± kv about Fermi momentum kp] this linearization 
is valid as long as Aq/ Ep ^ 1. From Fig. [T]it is evident that over a certain (higher) range 
of velocities v there are two solutions for A, these solutions differing from those obtained by 
Bagwell j^, by whom multiple solutions were not found. This type of multiplicity feature 
was first observed by Sarma ^ in the context of the exchange field in superconductors. In 
fact, the self-consistency conditions holding in the exchange-field case are identical to the 
ones holding here, provided one makes an identification between kpv and the exchange-field 
energy fiB h. The two situations appear to be essentially equivalent. 

The feature of double solutions that we have found for the order parameter, as shown 

n 

in Fig. [H is in contrast with previous findings [6j (see, in particular. Fig. 3 therein). We 
suspect that this discrepancy results from the use of an iterative scheme for solving the self- 
consistency equation numerically; when multiple solutions exist, such scheme yields results 
that depend sensitively on the initial conditions. We have instead used the bisection method, 
which locates all possible solutions. 
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V. TRANSITION TEMPERATURES: LINEAR INSTABILITY 



If one assumes that the transition between the superconducting and normal states is 
continuous, one can determine the transition temperature via the self-consistency condition 
taken in the limit of vanishing A. We note that, owing to the relation of the Fermi function 
to the Matsubara sum (see, e.g., Ref. [sl), Eq. (fTOj) is equivalent to 

"^"^^^^ {iuJn-kv-Ek){iuJn-kv + Eky ^^^^ 

where uJn = 27iT(n + (1/2)) with n = 0, ±1, ±2, . . . , are the Matsubara frequencies. We 
denote by the value of the temperature that solves Eq. 0151) in the limit A — 0. In 
this limit, imposing a cutoff ud on the Mastubara frequencies, and integrating out the 
momentum k, we arrive at the condition 

l = «iV„2.r.5:Re(— (16) 



OJn>0 



in which we have redefined Un to mean Un = 2TiTc{n + (1/2)). By exchanging g for the 
zero-flow critical temperature and performing the summation, Eq. f|T6|) becomes 

,„(|).,(i)_Ke„(i.|g). PT) 

where V'(^) is the di-gamma function [l^. 

The solutions of Eq. ( fT71) for Tc are shown in Fig. O which, like Fig. [1], exhibits a double- 
solution feature for a range of velocities. The solution at Tc = and kpv = Ao/2 corresponds 
to the branch A^ = 2kF v Aq — Aq. As we shall see below, this branch corresponds to an un- 
stable solution; a correct description of the vs. v phase diagram requires the consideration 
of free energies. 



VI. VELOCITY-TEMPERATURE PHASE DIAGRAM 



In the light of the results obtained in Sec. |Vl we ask the question: What is the true 
transition temperature? To determine this, one has to take into account the multiple solu- 
tions of the order-parameter self-consistency condition and compare their free energies (e.g., 
see Ref. ^). With this prescription we ask: Will our transition temperature Tc{v) have 
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the same functional form as the Tc{h) of Sarma 8] and of Maki and Tsuneto (SMT), in 
which h is the magnetic exchange field? 

In the case considered by SMT, a comparison is made between the free energies of a spin- 
unpolarized superconducting state and a normal state that is partially spin-polarized, both 
states being subject to an exchange field. In the present case the comparison is between a 
flowing superconducting state and a stationary normal state. Viewed from a reference frame 
that flows with the superconducting state, the normal state flows and can be regarded as 
the analog of SMT's polarized state. Furthermore, our superconducting state is the analog 
of SMT's superconducting state. This perspective suggests that the correct procedure for 
our purposes is to compare the free energies of the flowing superconducting state and the 
stationary normal state, as was first done by Bagwell 6|. 

The physical reason for the equivalence between the exchange-field case considered by 
Sarma and the superflow case considered here is as follows. In the case of Sarma, due to 
the exchange field, the up and down electrons in a Cooper pair have an energy difference of 
fisB. In the case of flow, the electron pairing is between states (k + q) t and {—k + q) |, and 
which have an energy difference of 2kFq/'rn near the Fermi surface. In the former case, the 
normal state can be polarized, and thus can reduce its energy by an amount proportional 
to h"^ . In the latter case, the stationary normal state has an energy that is lower than the 
flowing normal state by an amount proportional to q^. Hence, the two scenarios appear to 
be equivalent. 

To understand the root of our state-selection procedure, we begin with an analogy. Con- 
sider the liquid and crystalline states of a particular material. Following Callen 171], we 
observe that the information sufficient to specify uniquely an equilibrium state is greater for 
the crystalline state than for the liquid state, owing to the spontaneously broken translational 
symmetry of the crystal and the resulting low-energy Goldstone field, i.e., the displacement 
field. Thus, when assessing the relative thermodynamic stability of a liquid and a crystalline 
state, one must specify a displacement field for the crystal but one must not (and indeed 
cannot) do so for the liquid state. It is meaningful to speak of an equilibrium crystalline 
state that supports a specified static shear stress, but not so for a liquid, as stress would 
induce a nonequilibrium dissipative steady state of shear flow, in which entropy would con- 
stantly be being produced. In summary, there is a unique equilibrium state for the liquid, 
but there is a family of equilibrium states for the crystal, differing in their displacement field 



10 



(or strain field, or its thermodynamic conjugate, the stress field). One is entitled to consider 
the relative thermodynamic stability of a liquid and any one of this family of crystalline 
states. Returning now to the case of superconductivity, the precise analogy is on the one 
hand between the crystal and the superconductor, and on the other hand between the liquid 
and the normal metal. 

One is thus entitled to consider the relative thermodynamic stability of a stationary 
normal state and any one of the family of the fiowing superconducting states. That there 
is a family of flowing superconducting states originates in the spontaneously broken gauge 
symmetry of the superconducting state, which amounts to a controllable thermodynamic 
field, (i.e., the phase field), which is the analog of the displacement field of the crystal. 
(The analogy is: phase 4^ displacement ; velocity 4^ strain; current density 4^ stress.) By 
contrast, the normal state is unique. 

To obtain the phase diagram, one should then compare the free energy Fs(q,fi) of a 
fiowing superconducting state to the free energy F„(0,yu) of a stationary (i.e., not fiowing) 
normal state 6|]. The energy difference between the fiowing and stationary normal states 
plays the role of the paramagnetic energy. We shall see that the true equilibrium transition 
temperature then has exactly the functional form obtained by Sarma |8| and by Maki and 



Tsuneto 
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The free energy (Q of the superconducting state at a fixed chemical potential fi is given 
explicitly by 
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-E 



\2m 2m 



2 cosh /3 



V2m 2m 



A2 + 2 cosh. (3kv 
(18) 

The stationary, normal state at chemical potential — {q^ /2m) thus has free energy 



2m) 



\2m 



2cosh/3f— + i- 

\2m 2m 



• (19) 



Fn[0,fi 

zm,/ 

k k 

Our goal is to find the difference between two energies at the same chemical potential, which 

is given by 

^2 



F,(g,/i)-F40,/i) 



2 

F,(g,/i)-F40,/i-|-) 



F„(0,/i)-F„(0,/i-|-) 
V 2mJ 



(20) 



where the term in the second pair of square brackets equals —Nq^ /2m (i.e., the analog of the 
paramagnetic energy), with being the total number of electrons. Hence, from Eqs. (fT9|) 
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FIG. 3: Free-energy difference Kj(g = mv,iJL) — Fn{Q,ii) = Nq^ /2m — 5F (in units of A^^q^o) 
various temperatures T = (0, 0.2, 0.4, 0.85) (black dashed, blue dash-dotted, green dotted, black 
solid, respectively) vs. superfluid velocity v (in units of vi). When Fs{q,fj.) — F„(0,/i) is negative, 
the corresponding superconducting solution is stable. The red dot (D) indicates the transition 
point at Vc = vl/\/2 and = 0; see also the point D in Fig. HI The two branches for lower 
temperatures (T = 0, 0.2, 0.4) represent double solutions for the order parameter. The range of v 
for which a nonzero solution exists shrinks as the temperature increases, and the feature of double 
solutions disappears at sufficiently high temperatures; see also Fig. [TJ The inset shows a blowup 
in the range [0.35,0.7] for t;. 

and (120|) we arrive at the result that 



F,(g,/i)-F„(0,/i) 



iV^ + — -rVin 

2m g ^-^ 



iV^ + — -Ty In 
2m g ^-^ 



cosh /SE'fc + cosh fikv 



cosh /^Sfc + 1 



(21) 



What we should do to obtain the phase diagram is to compare the possible nonzero solu- 
tion(s) for the order parameter to determine which solution possesses the lower free energy 
Fs{q,fi) and whether this free energy is lower than that of a stationary (i.e., nonmoving) 
normal state, i.e., whether Fs{q, n) —Fn{0, n)<0. For convenience, let us make the definition 



SF{v, fi) = - F,{q, fi) + F„(0, /i) 

2m 



A2 

-— +TVln 
9 ^ 



cosh + cosh pkv 
cosh + 1 



(22) 



The solution for A that has the largest value of SF — Nq^ /2m, provided this difference is 
greater than zero, is the true equilibrium solution; otherwise the normal state will be the 
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FIG. 4: Temperature versus flow- velocity phase diagram, indicating normal (N) and superconduct- 
ing (S) regions. The phase boundary ABD indicates the transition between superconducting and 
normal equilibrium states, which is obtained using the principle described below Eq. (j2ip . Across 
segment AB of the phase boundary the superconducting-normal transition is continuous. Across 
segment BD the transition is discontinuous. As discussed in Sec. IVIH line segments BC and BE 
indicate linear stability limits for, respectively, the normal and superconducting metastable states. 
Tc and v are measured in units of (i.e., the zero-flow critical temperature) and vl, respectively. 

true equilibrium state. The kinetic flow term Nq^ /2m can be re-expressed (by relating the 
total electron number to the density of state) as A'^q Aq (v/vl)"^, where Nq is the density of 
states per spin, and we have, for convenience, set the total "volume" of the system to be 
unity. Note that from Eq. (!22|) we have at zero temperature 

6F{v, i^) = No^ + 2j2{kv-E, + ^)o{kv - Eu). (23) 

k ^ 

This expression also holds for higher-dimensional systems, provided we replace k hy v 
by V, and kv hj k ■ v. The effective condensation energy 6F for the solution A = Aq 
[see the discussion below Eq. (1131) ] is NoAl/2. Balancing 6F against the additional energy 
Nq^ /2m = Nq Aq {v /v^Y due to the flow, we obtain Vc/vl = 1/ \/2- This corresponds to the 
Clogston-Chandrasekhar limit in the exchange-field case 

In Fig. [3l we plot 5F versus the momentum q for the two superconducting solutions at 
several temperatures, as well as Nq^ /2m. The crossing of the curve representing Nq^ /2m 
by the curve representing the superconducting state signifies a discontinuous transition. By 
finding such transitions at various values of we arrive at the equilibrium phase boundary 
Tc(t;), shown as the curve ABD in Fig. HI 
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VII. SUPERCOOLING AND SUPERHEATING 



As discussed in Sec. |VT1 we have discussed the true equiUbrium phase transition curve 
ABD; see Fig. m But what is the physical meaning of the branch BC? It would be a limit 
of metastability if a normal state with flow were a true equilibrium state. Suppose that 
the system were in an equilibrium normal state with v > Vc (i.e., lying to the right of the 
equilibrium phase boundary ABD), and were then 'quenched' by the rapid change of T 
and/or v into the metastable region BCD. The system would remain in the normal state, 
but metastably so, until a fluctuation were to occur that would nucleate a droplet of the 
superconducting phase large enough to grow and complete the conversion of the state to 
the true equilibrium state, i.e., the superconducting state. The metastability limit (which is 
also known as a spinodal line) is formally obtained by solving for the temperature obeying 
the equation [c.f. Eqs. ([lO]), ([IS]) and ffT6|)] 

1 SF 

hm = 0- (24) 

However, as we argue in Sec. IVH a normal state with a constant flow is not a true equilibrium 
state; therefore, strictly speaking, the curve BC does not represent a meaningful supercooling 
curve. 

On the other hand, the branch BE is a true limit of metastability for the superconducting 
state. The system can be 'quenched' into a metastable superconducting state from an 
equilibrium superconducting state [i.e., from (T, t>) to the left of the equilibrium phase 
boundary] by rapidly changing (T, v) to a value in the region BED. The system will remain 
superconducting, but metastably so, until a normal-phase nucleation event carries it to 
the equilibrium (i.e., normal) state. In this case, the metastability limit is determined by 
simultaneously solving the equations 6F/6A = and 6'^F/6'^A = (but A ^ 0) numerically; 
see Fig. HI 

We note that a corresponding diagram in the case of exchange-field effect in supercon- 
ductors was first obtained by Maki and Tsuneto j^. To answer the question of how long 
it would take for system in a metastable state to find its equilibrium state is a kinetic one 
that would require further investigation. 
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FIG. 5: Supercurrent Iq (in units of Iq) vs. superfluid velocity v (in unit of vl) for various tem- 
peratures. From top to bottom: T = (0.1, 0.25, 0.4, 0.5, 0.556, 0.75, 0.9) r°. The curves terminate 
at the critical velocities Vc{T) appropriate to these temperatures. The maximum supercurrent for 
a particular curve determines the value of the critical current at that temperature. 

VIII. SUPERCURRENT 

In this section we examine the dependence of the supercurrent on the flow velocity, and 
thereby determine the critical current, i.e., the maximum equilibrium supercurrent that the 
system can sustain at various temperatures. In general, the charge current Iq carried by 
the system is defined via 

^« = ^ 5Z (4+5,^ • (25) 

fc,cr 

As E.,.( c^_i_g^o-Cfc+i3,o-) is the total number of electrons , we can re-write Eq. ([25]) as 

^ = ^ AT ^; + e ^ A (4^^^^ , (26) 

where v = q/m. By using Eqs. (jlj) and ([8]), the second term can be simplified to 
(2e/m) J2k ^ fT-;k, and thus the expression for Iq becomes 

Iq{T,v) = -\^\Nmv + 2y2kuA, (27) 

k 

where the electron charge is e = — |e| and the temperature T is implicit in the Fermi function 
see Eq. (E]). 

We shall compute the current in the superconducting state as a function of T and f , and 
then, by choosing the value Vm{T) of v in the range [0,fc(^)] that maximizes Iq{T,v) at 
fixed T, we shall obtain the critical current, denoted by Jc(T)[ = Jq(T, t>m(T))] , at that 
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FIG. 6: Critical current Ic (in units of Iq) vs. temperature T (in units of T^). This is also current 
vs. temperature phase diagram. 

temperature. Equation ( !27l) . applied at T = 0, tells us that the zero-temperature current is 
given by 



where Iq = — A^|e|Ao/A;^ = — 2|e|Ao/vr. Although, formally, Iq{T = 0,v) achieves its maxi- 
mum value (viz., Iq) at kpv = Aq, this maximum value is unattainable, because the super- 
conducting state gives way to the normal state via a pre-emptive transition at kpv = Ao/\/2 
(unless the system falls out of equilibrium and remains metabstably in the superconducting 
state). Therefore, the maximum equilibrium critical current at zero temperature occurs 
at V = vl/V^, and has the value Iq/^/2. At v = vl/V^ the superconducting state first 
becomes unstable, so the critical current that is attainable in practice is slightly less than 
Iq/\/2. This is unlike the behavior in three dimensions, as well as the higher of the nonzero 
temperatures; for these situations the maximum current is achieved at a value of superfiuid 
velocity for which the superconducting state is still globally stable. 

In Fig. [5] we show the equilibrium supercurrent Iq vs. the superfiuid velocity v for 
various values of the temperature. The value of f = Vm that corresponds to a maximum in 
the current (i.e., the critical current Ic) is not necessarily v^. (It is for lower temperatures, 
but not for higher ones.) The temperature dependence of the critical current Ic{T) is shown 
in Fig. El 



lQiT = 0,v) 



^ \Nmv - 2No kpJklv^ - AlQ{kFV - Aq) 



(28) 




(29) 
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FIG. 7: Superfluid density ns{T) (in units of 2ekp/TT) vs. temperature T (in units of T^) 

A related quantity of interest is the superfluid density n^, defined as 

djs{v) 



dv 



v=0 



(30) 



i.e., the response of the current density js (whose expression is identical to Iq because we 
have set the system 'volume' to unity) to an infinitesimal change in the 'driving' velocity 
V. The temperature dependence of Ug is shown in Fig. [71 The vanishing of near Tc is, 
as expected, linear, and the departure of Ug from its zero-temperature value has, at low 
temperatures, a thermally activated form. 

IX. EFFECT OF DISORDER 



So far, we have discussed clean systems, i.e., systems not disordered by any scattering 
electrons by impurities. What effect will disorder have? In particular, would the disorder 
change the nature of the superconducting-to-normal transition? In this section we briefly 
discuss the latter issue. In the case of disorder, the Green function technique is better 
suited for dealing with disorder than is trying to diagonalize the Hamiltonian with inclusion 
of an arbitrary configuration of the impurities. An equivalent calculation in the context 



of the exciange field was done by Maki and Tsuneto 



9| and, more recently, by Wei and 



Goldbart ISj in the context of Little-Parks effect in small rings. Here, in contrast with 
Ref. 1^, we are not concerned with the finiteness of the system size, and we simply quote 
the relevant results for the equation obeyed by the critical temperature, i.e., Eq. (B39) of 
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FIG. 8: The impact of elastic scattering on the flow- velocity-dependent normal-to-superconducting 
transition temperature solutions obtained via the self-consistency, Eq. (I31|) . and v are in units of 
and vl, respectively. From left to right: K/tq Aq = (0, 1, 1/0.546, 5, 10), i.e., cleaner to dirtier. 



Ref. 181: 



In 



-a + \/a^ - 



1 + a + - 



a + \/a^ -x^ , / 1 + a - \/ - x^ 



(31) 



where, in the present case, a = l/AnroTdv) and x = krv/nTc^v), the parameter tq is the 
elastic mean-free time, and ip{x) is the di-gamma function 16[ |. 

In the clean limit (i.e., tqAq ^ 1), Eq. ( ISTll reduces to Eq. ( ITTl) . for which multiple 
solutions for Tc{v) exist. For strong disorder (tq Aq <C 1), Eq. (I3T!) reduces to 



In 



TJv] 



i' 



(32) 



for which no multiple solutions for T^iv) exist. This suggests that the assumption of a 
vanishing order parameter in the search for the critical temperature is valid and, hence, 
that the transition is continuous. Moreover, from Eq. (132|) . in the strong-disorder limit the 
critical velocity is determined via 



4ro Aq. 



(33) 



Numerically, we have found that for tq Aq < 0.55 multiple solutions do not occur. This, then, 
gives the threshold for the disorder strength that divides the discontinuous and continuous 
transition regimes. In Fig. [8] we show the solutions of Eq. ( 13T]) for various values of tq Aq. 
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In the strong- disorder limit we obtain from Eq. (l33l) that Vc/ f l = \/ '^^o / 4/e (where le is the 
elastic mean-free path and is zero-temperature coherence length), in particular that the 
critical velocity exceeds v^. But will the critical current also exceed its clean-limit value? 



According to Bardeen {4], the superfluid density in the strong-disorder limit is reduced by 
a factor of roughly le/^o- An order-of-magnitude estimate gives for the critical current the 
product of the superfluid density rig ~ {h/^o)n^ (where denotes the superfluid density 
in the clean limit) and the critical velocity Vc = ^^/tt^oJ^vl. This, in turn, gives that the 



critical supercurrent is reduced by roughly a factor \/le/^o, which is much less than unity 
in the strong-disorder limit. Therefore, even though the critical velocity is increased by 
disorder, the critical current is reduced by it. 

We now digress to discuss the scenario of a closed superconducting loop threaded by 



magnetic flux, which is pertinent to the exp eriment performed by Liu et al 



situation considered theoretically in Ref. 



• 3 



and the 



18l ] . There, the boosted velocity is induced by the 
magnetic flux. The measurement of Tc vs. flux = $/$o (with $0 = hc/2e) presented in 
Fig. 4 of Ref. [isl shows non-$o-periodic behavior, with the first dome (near = 0) being 
higher than the second dome (near 0=1). Naturally, the data explore only a limited range 
of flux, so it is not possible to determine with certainty whether higher flux values would 
yield dome heights that oscillate with flux or decay. Oscillations would be a consequence of 
flux through the hole; decay would be a consequence of flux through the sample, which would 
cause pair breaking. To ascertain which of these two effects (oscillation or decay) is likely 
to dominate in the range of fluxes probed by the experiment, it is useful to compare several 
lengthscales. First, the sample radius R ~ 75 nm is smaller than the coherence length C,o, by 
a factor of roughly 0.06. Following Ref. |l^, this suggests the smallness of the radius as the 
origin of the dome-height mismatch. On the other hand, a more complete characterization 
of the system would take into account the elastic mean-free path l^, which for the sample 
studied in Ref. [l^ is shorter than the radius, by a factor of approximately 7.8. This puts 
the sample in the dirty regime, and as a result, the impact of the smallness of the radius is 



suppressed and tends to restore the $0 periodicity of Tc discussed in Ref. [l8|] and given by 

,„(ZM)..,,,.,(1)-,(1.£^||M), (3, 

where F ^ 1.76 and Xm(0) = minmez |0 — m\/2. What else could cause the lack of $0 
periodicity? Because the magnetic flux is not confined inside the hole in the loop, but is also 
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FIG. 9: Critical temperature Tc{cl)) (in unit of zero-flux critical temperature T^) vs. magnetic 
flux [(f) = <I>/(/ic/2e)] for a superconducting ring. The result is obtained by using Eq. (I34p with the 
additional pair-breaking term given in Eq. (|35p . and a numerical factor of about 2.5 is the only 
fitt ing parameter. All other parameters are taken from those reported in the experiment by Liu et 



al. 



present throughout the sample, owing to its finite thickness, the orbital pair-breaking effect 
comes into play. To account for this effect, we use the scenario of the orbital pair-breaking 
effect of a parallel magnetic field applied to a film of finite thickness d (which is about 25 nm 
in the experiment), which leads us to add to the argument of the second term on the r.h.s. of 
Eq. the term 

1 De^B^d^ T^\,led\ ^3^^ 



see, e.g., Ref. (20[. (The term should be correct up to a numerical factor due to geometry.) 



As we see from Fig. [9l the resulting dependence of the critical temperature on the flux 



seems to reproduce, at least quahtatively, the features reported by Liu et al. [l9|. 



X. HIGHER DIMENSIONS 



We have seen that in the clean limit the discontinuous transition at flow velocity v = 
Ao/kpV^ occurs due to a competition between the condensation energy A^'q A^/2 and the 



flow kinetic energy Nq^/2m = NqAKv/vl)^ (for one dimension; see footnote [211]), where 
vl = Ao/kp- The principle that we employed in Sec. lVIl is to compare the free energies of the 
flowing superconducting states and of the stationary normal state. As further confirmation 
of the validity of this principle, we now apply it to higher- dimensional settings, and compare 
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FIG. 10: The energy difference Fs{q = mv, ji) — Fn{0, n) (in units of A^'o Aq/2) versus flow velocity 

V (in units of vl) in three spatial dimensions. When the difference is less than zero, the flowing 
superconducting state is stable. 

the results with the known one mentioned in the Introduction. 

In two dimensions, the solution of the self-consistency equation fill I) gives A = Aq for 

V < Vl, and A = for v > v^. This indicates a discontinuous transition at vf'^^ = vl. 
Indeed, this result is consistent with free-energy considerations: by using the condensation 
energy 5F = Nq Aq/2 and the flow kinetic energy Nq^ /2m = NqAq{v/vl)^ /2 (the difference 
from the one-dimensional case arising from the density of states; see footnote 2]|), we 
conclude that = vl. 

In three dimensions the solution of the self-consistency equation (ITT]) gives A = Aq for 

V < Vl, and for vl < v < (e/2) vl (see Ref. js]) gives the following implicit equation for A: 



1 



A 

kp V 



In 



1 + A 1 



A 

kp V 



(36) 



where we have ignored terms of order A/ujl> or smaller (with ud denoting the Debye fre- 
quency). Equation ( l36i) gives that A = at f = v[^^^ = {e/2)vL ~ 1.359fL. Is this 
linear-instability result consistent with free-energy considerations (i.e., the balancing of 
the two energies)? If one were to naively use the condensation energy Nq/S.'^/2 to bal- 
ance the flow energy Nq^/2m = Nq A'^{v/vlY /'^ ^5 would obtain v^. 



v/372 



1.225) i;l < vf^^l Does this mean that the superconducting state really becomes unstable 
1.22b Vl < v^^^l Such a result would be in conflict with the result, mentioned 



at V 



in Sec. HI that 



,(3D) 



(e/2)t;L? How is this apparent conflict resolved? For v < vl, the 



effective condensation energy 6F is indeed A''oAq/2. But for v > vl, SF 7^ A''oAq/2, so 
it would be incorrect to equate NqAI/2 and Nq^/2m in order to deduce v^- Instead, one 
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should use the expression for 5F appropriate to the range v > vl- In fact, an evaluation of 
Eq. fl23|) in three dimensions gives (for v > vl) 

3 



5F= f A^ + f(W 



A 

kpv 



+^„A^,M,-.«H5!_A,„A.,„V(iZS!Z5i±^, (37) 

where the dependence of A on is given by Eq. (!36l) . From Eqs. (123|) and (!37j) we have 
that for Vl < V < {e/2)vL, A > and 5F — Nq^/{2m) > 0. Only when v increases to 

V = {e/2)vL (i.e., when A = 0) do we have 6F — {Nq^/2m) = 0. Thus, the transition at 

V = {e/2)vL is continuous, as the order parameter vanishes continuously there, and, hence, 
the conflict is resolved. The dependence of the quantity Fs{q, yu) — F„(0, /i) = Nq"^/ {2m) — 6 F 
on flow velocity v is plotted in Fig. [TOl which clearly shows that the superconducting state 
is stable for flow velocities in the range < v < {e/2) v^. 

XI. CONCLUDING REMARKS 

We have revisited the issue of the critical velocity of an effectively one-dimensional super- 
conductor using mean-field theory. As we have discussed, this issue is equivalent to that of 
the threshold uniform exchange-field in a superconductor, first investigated by Sarma 8|] and 
by Maki and Tsuneto 91]. This equivalence is due to the correspondence of (i) the Zeeman 
frustration energy fisB between the two electrons of a Cooper pair (in the case of exchange 
field), and (ii) the kinetic frustration energy 2kFV (in the case of flow). In particular, we 
have found that at zero temperature the critical velocity is a factor of smaller than 
the Landau critical velocity, in contrast to a previous finding j^. Our result originates in 
a Clogston-Chandrasekhar-like discontinuous transition between the superconducting and 
the normal state. The physical reason for the discontinuity of the transition is that it results 
from a balance between the condensation and flow energies. (Recall that the bulk supercon- 
ducting phase transition in magnetic field results from a balance between the condensation 
(A'oAo/2) and magnetic field tt) energies, and that the transition is discontinuous.) 

The transition remains discontinuous for temperatures below approximately 0.56T°; above 
that, it is continuous. 
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We have also studied the issue of critical supercurrents, determined the current- 
temperature phase diagram, and examined metastability (and its limits) in the temperature 
versus flow-velocity phase diagram. As a test of our underlying principle, namely the com- 
parison of the free energies of flowing superconducting and stationary normal states, we have 
also examined the two- and three-dimensional cases, and have thus obtained results that are 
in agreement with previous findings. In addition, we have commented on the effects of 
electron scattering by impurities and, in particular, we have argued that strong disorder 
renders continuous the aforementioned discontinuous phase transition. Even though disor- 
der can increase the value of the critical velocity, the physically measurable quantity, i.e., 
the critical current, is still reduced by the disorder. 

Throughout this Paper, we have adopted a mean-field approach, thus ignoring the effects 
of fluctuations such as phase slips. Furthermore, we have limited ourselves to the analysis 
of equilibrium properties as well as issues of metastability. The intriguing issue of how 
a flowing superconducting state undergoes the transition to the stationary normal state, 
which essential dynamical processes take place, and on what timescale all remain as research 
directions for the future. 
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